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Abstract. We compare the one-loop renormalization group flow to strong coupling of the electronic inter- 
actions in the two-dimensional t-t'-Hubbard model with t' = — 0.3t for band fillings smaller and larger than 
half-filling. Using a numerical iV-patch scheme (N = 32 . . . 96) we show that in the electron-doped case 
with decreasing electron density there is a rapid transition from a d x 2_ y 2-wa,ve superconducting regime 
with small characteristic energy scale to an approximate nesting regime with strong antiferromagnetic ten- 
dencies and higher energy scales. This contrasts with the hole-doped side discussed recently which exhibits 
a broad parameter region where the renormalization group flow suggests a truncation of the Fermi surface 
at the saddle points. We compare the quasiparticle scattering rates obtained from the renormalization 
group calculation which further emphasize the differences between the two cases. 
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1 Introduction 

Viewed from a weak coupling perspective, the two-dimen- 
sional Hubbard model exhibits an interesting interplay be- 
tween different types of fluctuations. Particular attention 
has been devoted to the situation close to half band fill- 
ing, mainly because of its relevance to the high- 2^, super- 
conducting cuprates jlj. For these particle densities, be- 

, cause of the approximate nesting between opposite sides 
of the Fermi surface (FS) or due to the large density 
of states around the saddle points of the dispersion at 
(7r,0) and (0,7r), the scattering processes with momen- 
tum transfer « (tt, tt) will be strongly enhanced. Besides 

, the fact that these scattering processes drive antiferro- 
magnetic (AF) fluctuations, they also create a d x 2_ y 2- 
wave component in the pair scattering Q . One-loop renor- 
malization group (RG) techniques represent a powerful 
method to analyze the coupling between antiferromag- 
netic and superconducting fluctuations in an unbiased way 

(see [§§|4ili-0'[yJ and thc P resent paper). 

Apart from the competition between two different types 
of ground states, the hole-doped t-t' Hubbard model close 
to half filling may exhibit a richer variety of possible ground 
states especially when umklapp scattering between elec- 
trons becomes important. In M and it was argued 
that in the case of not-too-small values of the next-nearest 
neighbor hopping t' , e.g. t' = — 0.3i, when the FS is close 
to the saddle points at (ir, 0) and (0,7r), the RG flow to 
strong coupling between the k-space regions around the 
saddle points has strong similarities with the RG flow in 
the half-filled two-leg Hubbard ladder. In the latter sys- 
tem, the ground state is well understood and is an insu- 
lating spin liquid (ISL) fl^| . Although there is no reliable 
theory for the two-dimensional case, the similarity of the 



RG flows suggests that here as well the true strong cou- 
pling state of the k-space regions around the saddle points 
is an ISL with a truncated FS. The parameter region in 
which this interesting flow to strong coupling occurs was 
called the saddle point regime. 

In this paper we investigate the RG flow to strong 
coupling for band fillings larger than half-filling, corre- 
sponding to electron-doping the half-filled system. Again 
we choose t' = — 0.3t. The main observation is that in the 
electron-doped case the saddle point regime of the hole- 
doped case is absent. Instead of the mutually reinforcing 
flow between AF and d-wave pairing processes which was 
found in the saddle point regime, there is a clear separa- 
tion between the energy scales of the two channels. We do 
not find any signs for a Fermi surface truncation around 
the saddle points. Moreover the "hot" Brillouin zone (BZ) 
regions responsible for the leading flow are now located 
closer to the BZ diagonals. The different result is a direct 
consequence of the location of the Fermi surface, which in 
the electron-doped case crosses the umklapp surface (US, 
defined in Fig. [|) in the BZ diagonal, therefore giving rise 
to strong scattering between the FS segments on oppo- 
site sides connected by the wave vector (tt, tt). These AF 
processes generate only a weak attractive d x 2_ y 2 compo- 
nent in the pair scattering channel and do not drive the 
pairing processes at low scales as they do in the saddle 
point regime of the hole-doped case. Consequently, upon 
increasing the band filling such that the (ir,ir) processes 
get cut off, we obtain a rather abrupt transition from a 
flow to strong coupling dominated by umklapp and AF 
processes with high critical scale to a Kohn-Luttinger- 
like instability with predominant d x 2_ y 2 -Cooper processes 
at comparably low critical scale. Therefore the electron- 
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doped case considered here strongly resembles the hole- 
doped cases with small absolute values of t analyzed in 
H and and differs from the hole-doped case with 

r = — 0.3i discussed in flipf . 

As another piece of information we present results for 
the scattering rates of quasiparticles at the FS with a RG 
method described below. These calculations further illus- 
trate the differences between the hole- and electron-doped 
cases. According to model calculations by Ioffe and Mil- 
lis [^3| and also Hlubina most transport experiments 
of the optimally hole-doped cuprates can be properly de- 
scribed by assuming a pronounced anisotropy in the scat- 
tering rate. The anisotropy is such that the quasiparti- 
cles around the saddle points scatter strongly and their 
spectral weight becomes smeared out while the quasipar- 
ticles in the BZ diagonal are only subject to a weak Fermi- 
liquid-like scattering. Our RG calculations with frequency- 
independent coupling constants can certainly not give re- 
liable information about the frequency dependence of the 
selfenergy. On the other hand their k-space resolution 
is rather good. Applying the RG scheme we observe a 
pronounced angular dependence of ImZ^k^a; = 0) with 
maximal scattering for particles at the saddle points and 
weaker scattering in the BZ diagonals. In the electron- 
doped case, the anisotropy is weaker and the FS regions 
with the largest scattering rate are located in the BZ di- 
agonals, tied to the intersection of the FS with the US. 

Angular resolved photoemission (ARPES) allows one 
to measure the temperature dependence of Im S(u> = 
0, kp) directly from the width of the quasiparticle peaks. 
For optimally hole-doped Bi 2212, Valla et al. [|l5) found a 
linear-T dependence almost everywhere on the FS. For a 
comparison we calculate the quasiparticle scattering rates 
at high temperatures, where the RG flow does not diverge. 
Our results for the hole-doped case are qualitatively sim- 
ilar to the experimental data. In the electron-doped case 
the RG calculations yield a non-linear T-dependence of 
the scattering rates closer to conventional T 2 behavior. 



2 The method 

Here we describe the RG method and the scheme for the 
calculation of the quasiparticle scattering rates. Readers 
who are only interested in the results can proceed to the 
next section. 



2.1 RG flow of interactions and susceptibilities 

The interplay between the different fluctuations mentioned 
in the introduction can be appropriately analyzed within 
one-loop RG where one successively integrates out inter- 
mediate states according to their band energy e^. For 
the two-dimensional Hubbard model on the square lattice 
with nearest neighbor hopping t and next-nearest neigh- 
bor hopping t' , 

ek = — 2i (cos k x + cos k y ) — At' cos k x cos k y — fi , 




x=>cx 



Fig. 1. The one- loop RG equations for two - and four-point 
vertex function. The dot symbolizes the derivative with respect 
to the energy scale A. 

where \x is the chemical potential and the lattice constant 
is set to unity. 

The RG scheme for 1PI irreducible vertex functions 
which we use is explained and discussed in detail in Ref. 
[ [l6| . Here we give simple graphical explanation of the RG 
flow for the two- and four-point vertex functions. The flow 
of higher order irreducible vertex functions is not taken 
into account. 

The differential equations for the two-point vertex (yield- 
ing the scale-dependent selfenergy -£Vi(k, iu>)) and the four- 
point vertex function, which describe their RG flow with 
decreasing energy scale A, are described graphically in Fig. 
^. At the internal vertices we have four-point vertex func- 
tions at scale A connecting the incoming and outgoing 
lines with the internal ones. Crystal momenta and Mat- 
subara frequencies are conserved at the vertices. One of 
the two internal lines corresponds to a so-called single 
scale propagator Sa which is only non-zero for modes at 
scale A while the other line denotes a full Green's function 
G a with low energy modes below scale A cut out through 
a cutoff function C(A, |ek|)- Throughout this paper we ne- 
glect possible self-energy corrections, therefore we have 

C(A,|e k |) d A C(A, \e k \) 
Gyi(k, iuj) — and «5yi(k, iuj) — . 

llx) — £k IUJ — £k 

The scale-dependent spin-rotation invariant four-point 
vertex function /^(ki, s\, k2, s 2 , k3, S3, k4, S4) can be ex- 
pressed in terms of a coupling function V^(ki, k 2 , k 3 ) [|l0| 
p"6| which is determined for the case s% = S3 and S2 = 
S4. The Si denote the ^-components of the spins of in- 
coming (si and S2) and outgoing (53 and S4) particles. 
VA(ki,k 2 ,k 3 ) is represented graphically in Fig. g. This 
means that the electron-electron interaction at energy scale 
A is a function of the two incoming wave vectors ki and k 2 
and one outgoing wave vector k3. For the bare Hubbard in- 
teraction at starting scale A rj At, we set Va (ki , k 2 , k 3 ) = 
U . For all calculations in this paper we choose U = 3t. The 
flow of V}i(ki, k 2 , k.3) with decreasing the energy scale A 
is then given by all types of one-loop diagrams including 
particle-particle and particle-hole diagrams with all pos- 
sible momentum transfers (see Fig. ||). 

For the numerical integration of these coupled equa- 
tions we use a phase space discretization following Zanchi 
and Schulz ||. The BZ is divided up into nN patches cen- 
tered around N lines. Each line starts from the origin in 
a certain angular direction 9(k) and from the (±7r, ±7r)- 
points so that the lines meet at the umklapp surface. The 
phase space segments around these lines are then further 
split into n patches. Here we took n — 3 where one patch 
per line is centered around the FS with e.g. |e(k)| < OAt, 
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k2,s k^,s 



Fig. 2. The coupling function VX(ki, k2, IC3). ki,s (k.2,s') 
specify the first (second) incoming particle, k3,s (k4,s') be- 
long to the first (second) outgoing particle. 




(e) 



Fig. 3. The contributions to the right-hand side of the 
RG equation for the coupling function Vyi(ki, k2, ks). (a) the 
particle-particle term (b) the crossed particle-hole term (c) the 
direct particle-hole terms; the first of these three graphs gets 
a factor —2 because of the fermion loop. 

while the other two patches cover phase space regions at 
higher positive or negative band energy. 

Then the coupling function is discretized as follows: 
we approximate Vyi(kx, k 2 , k 3 ) by a constant for all wave 
vectors in the same patches and calculate the RG flow 
for the subset of interaction vertices with one wavevector 
representative for each patch. In the low energy patches 
around the FS we choose to take these wave vectors as the 
crossing points of the N lines with the Fermi surface (FS), 
the wave vectors for the remaining patches away from the 
FS correspond to a higher band energy, e.g. |e(k)| = 0.8t 
(see Fig. ^). The phase space space integrations are per- 
formed as sums over the patches and integrations over the 
radial direction along 3 or 5 lines inside each patch. Most 
calculations were done using a 32x3 system, i.e. N = 32 
and n = 3. Two typical Fermi surfaces with TV = 32 points 
are shown in Fig. pi 

Along with the flow of the interactions we calculate 
several static susceptibilities by coupling external fields of 
appropriate form to the electrons. During the flow these 
external coupling are renormalizcd through one-loop ver- 
tex corrections, as described in Ref. jn^. The pairing and 
AF susceptibilities can be calculated directly with the RG 
flow of the coupling function. The susceptibilities and cou- 
plings to uniform external charge and spin fields require 
special care as they only involve particle-hole excitations 
in a shell around the FS with width of the temperature. 
Before the cutoff A reaches that scale the interactions have 
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Fig. 4. The Brillouin zone, Fermi and umklapp surfaces and 
the lines in the patch centers for N = 32. The solid dots in 
the patches denote the wave vectors for which the four-point 
vertex function is calculated. For all wave vectors inside the 
same patch, VSi(ki, is approximated by a constant. The 

solid line denotes the non-interacting FS for t' = — 0.3i and 
fi = —t, which is in the saddle point regime of hole-doped 
case. The curved dashed line is a typical FS for the electron- 
doped case. The straight dashed lines connecting the (±7T, 0)- 
and (0, ±7r)-points denote the umklapp surface. The arrows 
symbolize umklapp processes between the saddle point and 
the BZ diagonals, respectively. 




Fig. 5. The Fermi surfaces and N = 32 points ki?(fci) for 
which the coupling function Vk(fci, &2, £3) is calculated. Left 
plot: fi = -1.05* (hole-doped, (n) « 0.83). Right plot: fi = 
(electron-doped, (n) ~ 1.2 per site). 



typically grown larger than the bandwidth such that the 
one-loop flow is not reliable any longer. In order to get a 
qualitative result for these cases we stop the flow of the 
interactions at a certain scale Af Teeze and treat them as 
constants in the continuation of the flow of the suscep- 
tibilities down to A = where we obtain the main con- 
tributions to the uniform susceptibilities. This is basically 
equivalent to stopping the RG flow at A[ rcczc and calculat- 
ing these susceptibilities for the effective theory below the 
cutoff /lfroczc within RPA (random phase approximation) 
using the renormalized interactions given by the RG down 
to ylfreczo- A similar way to obtain the uniform susceptibil- 
ities from the flow of the Landau function has been used 
in Ref. §. 
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2.2 Approximate calculation of the scattering rates 

Here we describe an approximation scheme for the calcula- 
tion of the imaginary part of the selfenergy from the flow of 
coupling function. The technical problem occurring in the 
RG formalism presented above is that we have neglected 
all frequency dependencies of the four-point vertex, hence 
the four-point vertex in the one-loop contribution to S at 
scale A in Fig. |l| is a real quantity and the diagram does 
not give an imaginary part. This deficiency can be repaired 
to a large extent by replacing the approximate frequency- 
independent vertex with the solution of the one-loop RG 
equation for the four-point vertex in its integral form. This 
effectively yields a two-loop term corresponding to a two- 
particle-one-hole intermediate state which, after analyti- 
cal continuation to real frequencies, gives a non- vanishing 
imaginary part. Moreover these contributions to Im Ua 
are the leading ones, because all neglected contributions 
from further insertions of one-loop diagrams at the ver- 
tices have intermediate states with a number of interme- 
diate particles larger than 3 and can therefore be expected 
to give only small contributions due to phase space restric- 
tions. If we ignored the flow of the coupling constants and 
integrated the expression thus obtained down to A = 0, 
we would retrieve the bare two-loop selfenergy. In this pa- 
per we neglect the feedback of the selfenergy on the flow 
of the coupling function. In principal this effect is con- 
tained in the RG equations but keeping it would increase 
the complexity of the calculations considerably. 

Formally we can calculate the selfenergy Sa=o()^f, *w) 
at the FS for Matsubara frequency ito by integrating the 
RG differential equation for the scale-dependent selfenergy 
depicted in Fig. p] as follows: 




[-2K 1 (k J ,,k',k F ) + V^(k F ,k',k')]^(k',z W ') 

For Va we now substitute the differential equation for the 
four-point vertex, integrated from the starting scale Aq 
down to A. This corresponds to inserting the one-loop 
diagram on the right side of Fig. [I] into the self-energy 
graph on the left side. Then we get a second integral over 
dA'; 

nO pA 

E A =o(kF,iu) = dA dA'Y,VA>GA>SA>S A V A > (1) 

JAo JA 

Here the sum is over all internal wavevectors, frequencies 
and spin indices of the different possible diagrams. For 
simplicity we have suppressed the arguments of the four- 
point vertices and propagators in this expression. Obvi- 
ously, repeating the insertion of scale-integrated one-loop 
terms for the vertices yields a perturbation expansion of 
the self-energy with arbitrary numbers of intermediate 
states. Here we restrict the analysis to the contributions 
from the two-loop term (|l|), which contains the Landau- 
Fermi liquid selfenergy and the deviations from it through 
the flow to strong coupling of the coupling function Va- 



(a) (b) (o) 




Fig. 6. The different diagrams for the two-loop selfenergy. Di- 
agrams (a) and (b) contribute to the imaginary part of the 
selfenergy, contributions of type (c) are real for external fre- 
quency uj + iS. There are 3 other diagrams of type (c) with 
different orientation of the vertices which are not shown here. 

In the integrand in (0) , the two four-point vertices Va' , 
one single scale propagator Sa> and one full propagator 
Ga> depend on A'. The only yl-dependence is in the differ- 
entiated cutoff- function C(A) of the original single-scale 
propagator Sa- Therefore (Q) has the structure 

n0 pA 

E A =oQtF,iu) = ^2 dAC(A) / dA' R(A') . 

JA JAo 

The double scale integral is numerically expensive, but we 
can circumvent it by a partial integration, resulting in 

r° 

£a=oQlf, i") = J2 dA - C ( A )] R ( A ) ■ ( 2 ) 

JA 

The remainder R(A) contains one single scale propagator 
at A and a propagator for modes above A. The difference 
of cutoff functions in the rectangular brackets means that 
the third of the three internal lines has its support on 
modes below the cutoff. The form (||) only contains a sin- 
gle yl-integral and the vertex functions at A. Therefore it 
can be integrated along with the flow of the four-point ver- 
tex. Note however that for consistency with (|l]) we should 
always integrate out the full scale range. 

Diagrammatically the insertion of the one-loop flow 
into the one-loop diagram for the selfenergy yields three 
topologically different diagrams which are shown in Fig. |(]. 
Each of these three diagrams gives 6 contributions: one of 
the internal line contains modes above A, another modes 
at A, and the third below A. 

Since we are interested in the quasiparticle scattering 
rates, we perform an analytical continuation onto the real 
frequency axis: the integrand of the two diagrams (a) and 
(b) in Fig. || contains the factor 

1 

iuj - ei - e 2 + £3 

which upon replacing ilu — > u> + iS with to = yields an 
imaginary part oc (5(e 1 + e 2 ~ £3)- The 5-function is smeared 
with a small width for the numerical treatment. 

One can show that for scale-independent vertices, e.g. 
all V/t(ki, k 2 , ka) = U, this scheme is equivalent to the cal- 
culation of the bare two-loop diagram. Thus for ordinary 
cases where the flow does not diverge we should obtain 
two-dimensional Fermi- liquid results. In our case the flow 
goes to strong coupling, and this leaves two possibilities 
for the analysis: either we choose a high temperature such 
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that the couplings do not become too large, then we can 
apply the above scheme down to zero scale and obtain an 
estimate for the imaginary part of the selfenergy above 
the strong coupling phase. The other option is, if we want 
to analyze the situation in the strong coupling regime, to 
stop the flow of the couplings at some scale Af ICCZC and 
integrate the flow down to zero scale with fixed couplings. 
We can then vary Af rcczc and thus obtain the change in 
the imaginary part of the selfenergy due to the flow of the 
couplings. This will tell us on which FS parts the renor- 
malized interactions weaken the quasiparticles most and 
we will e.g. find that in the hole-doped case the quasi- 
particles in the BZ diagonals are not affected by strong 
scattering. The results of both types of analysis are de- 
scribed in Section |[ 



3 Comparison of the flows to strong coupling 

First we compare the RG flow of the interactions. For the 
chosen parameters we generally find a flow to strong cou- 
pling, i.e. at sufficiently low scales and temperatures some 
components of the coupling function VA(ki, fe, fe) take 
absolute values larger than the band energy. We analyze 
the flow to strong coupling at scale Aw where the coupling 
functions just exceed the order of the bare bandwidth, i.e. 
at VA w ,max. ~ 8 — 12t. At these scales and for typical tem- 
peratures, the coupling functions have already developed 
a pronounced k-space structure and the dominant interac- 
tion terms at this scale will certainly be important for the 
strong coupling state. On the other hand the FS shift due 
to one-loop selfenergy corrections is still small |T(J such 
that it does not qualitatively change the flow above Aw- 
Similarly the scattering rate for particles at the FS (see be- 
low) remains smaller than Aw- Note however that at Aw 
and for initial interaction U = 3t the flow has not reached 
an asymptotic form and different classes of coupling con- 
stants would evolve differently if we continued the flow 
below Aw (where our method breaks down). Therefore 
the analysis of the flow to strong coupling remains qual- 
itative and does not provide definitive conclusions about 
the true strong coupling state. 



3.1 Hole-doped model 

The hole-doped case with t' = — 0.3i was extensively dis- 
cussed in Ref. [ ^0[ . In order to simplify the comparison 
with results for the electron-doped case we repeat the 
main observations briefly. We found that in the density 
region where the FS crosses the US close to the saddle 
points the dominant processes which grow fastest upon in- 
tegrating out higher energy modes are given by mutually 
reinforcing umklapp and Cooper processes between the 
saddle point regions. This leads to a flow to strong cou- 
pling at relatively high critical scales. The scattering pro- 
cesses involving quasiparticles in the BZ diagonals grow 
as well but they are merely driven by the couplings be- 
tween the saddle points and diverge less rapidly. These 
results are shown in Fig. ^ and Fig. |^. Starting the flow 



V/t 




8 16 24 



Fig. 7. Snapshot of the RG flow when the couplings have 
reached values larger than the bandwidth for the hole-doped 
case in the saddle point regime, t' — — 0.3t and fi = — t. The col- 
orbar denotes the value of the couplings V(ki, k2, ka)/t, where 
kF(fci) and kf(fc2) are the two incoming wavevectors on the 
FS and the outgoing wavevector kF(fo) is fixed at point 1 or 4. 
The 32 points are numbered according to their position around 
the FS, points 1,8,9,16 etc. are closest to the saddle points. 



with purely repulsive interaction, the Cooper channel will 
initially always decrease the values of the Cooper scat- 
terings. However if the FS is close to the saddle points 
the particle-hole channel with the rapidly growing umk- 
lapp processes also affects the large angle pair scattering 
processes which involve momentum transfer (7r,7r). In the 
saddle point regime of the hole-doped case this effect is 
stronger than the suppression through the Cooper chan- 
nel and the large angle pair scattering processes grow upon 
reducing the energy scale which in turn enhances the flow 
of the small angle pair scattering to negative values. Sim- 
ilarly it can be also seen that the growth of the d-wave 
pair scattering enhances the flow of the umklapp (tt, it) 
processes. In the saddle point regime with the FS close 
to the (0,7r) and (tt, 0)-points this mutual reinforcement 
between <i-wave and AF processes is very effective as it 
only involves low-energy quasiparticles. Therefore d-wave 
and AF processes diverge together at a single high energy 
scale. 



3.2 Flow in the electron-doped case 

For the electron-doped side we again take t' = —0.3t. 
Here, unlike the hole-doped case, smaller absolute val- 
ues of t' yield a similar picture, only the critical scales 
become larger due to the improved nesting. Half-filling 
corresponds to fi = — 0.65i. Upon increasing fi and (n), 
the FS crosses the US close to the BZ diagonals until 
/i = Ot and then loses contact to the US if we further 
increase the particle number. At these densities the flow 
to strong coupling changes drastically therefore we con- 
centrate on this filling range. In Fig. |lC] we show the flow 
to strong coupling for two different chemical potentials: in 
the upper plots [i = —OAt and the FS intersects the umk- 
lapp surface close to the BZ diagonals. Here the dominant 
processes which become large at comparably high energy 
scales are umklapp processes between the BZ diagonals 
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Fig. 8. Flow of certain coupling functions and susceptibilities 
in the saddle point regime of the hole-doped case with t' — 
— 0.3t and fi = — t. The dashed lines in the left plot show the 
flow of Cooper couplings, the solid lines umklapp processes 
between the saddle points. The dotted line denotes the largest 
coupling in the BZ diagonal. In the right plot the solid line 
denotes the a!- wave pairing susceptibility Xd-wave, the dashed- 
dotted line the AF susceptibility Xs{t,iv). 



involving the momentum transfer (it, tt) (see solid lines in 
Fig. and the sharp features in Fig. g). 

The Cooper processes between particles at the FS have 
dominant d x i_ y i -symmetry but remain of the order of the 
bandwidth even when we integrate the flow far out of 
the perturbative range. Thus the coupling to the strongly 
growing AF processes is only weak, as expected from the 
location of the FS. Consequently the characteristic en- 
ergy scales where AF and d-wave channels start to grow 
are very different, in strong contrast to the saddle point 
regime of the hole-doped case. Nonetheless we can still 
find signs of a coupling between d-wave and AF channel: 
if we consider Cooper processes for particles at the saddle 
points below the FS we find a much stronger growth to- 
wards lower scales (see dashed-dotted lines in Fig. [l(]). A 
large angle Cooper pair scattering between two inequiv- 
alent saddle points involves momentum transfer ss (n, tt) 
and is therefore enhanced by the AF processes. But un- 
like the saddle point regime in the hole-doped case, the 
fc-space regions where the coupling occurs are away from 
the FS and therefore there is only a weak influence of this 
mechanism on the true low-energy excitations. One should 
also note that in more realistic calculations including self- 
energy effects the coupling through off-FS processes will 
be further reduced by the short lifetimes of the quasipar- 
ticles in these regions. 

Increasing the band filling such that the FS loses its 
intersection with the US, the umklapp and other (tt, tt)- 
processes in the diagonals get cut off and saturate at low 
scales (see solid lines in the lower left plot of Fig. |l0|) . Thus 
the ci-wave Cooper processes can finally diverge at very 
low scales (at least two orders of magnitude lower than in 
the previous case). The Cooper pair scattering is shown 
in in the inset of Fig. [l0| Again the off-FS d-wave Cooper 
processes closer to the saddle points benefit more from the 
AF scattering processes and grow more strongly than the 
Cooper processes for pairs at the FS, but this difference 
becomes less with increasing particle density away from 






















































































































































































































































































Fig. 9. Snapshot of the RG flow when the couplings have 
reached values larger than the bandwidth for the electron- 
doped case, t' — — 0.3t, T = 0. For the upper two plots fi = 
-0.lt ((n) » 1.18), and the FS intersects the US. The left plots 
show the dependence of the scattering vertex Ui(fci, &2, &3 = 1) 
with the first outgoing wavevector fixed at point 1 on the in- 
coming wavevectors kf(fci) and on kj^fe) where ki labels the 
positions of the points around the FS (see Fig. |^). In the upper 
right plot, ks — 4, i.e. the hrst outgoing wavevector is near the 
BZ diagonal. The sharp vertical repulsive feature at fc2 = 21 is 
due to nesting with the wavevector kF(fe) — kF(fcs) ~ (n,n). 
For the lower plots, fi = 0.04f ((n) « 1.22) and the FS lies 
outside the US. The left plot again corresponds to kz = 1 and 
the right plot is for k$ — 4. The diagonal features belong to 
Cooper processes with kf(fci) + k^fe) — 0. The sign struc- 
ture of the Cooper pair scattering reveals that the dominant 
pairing symmetry is d x 2_ y 2, e.g. Va(1, 17, 1) is attractive while 
V/i(l, 17,9) is repulsive. 



half-filling. Apart from this side-effect we emphasize that 
in this low-scale e?-wave regime only Cooper processes di- 
verge. Umklapp and and nesting processes are cut off and 
do not develop singularities. 



4 Flow of the susceptibilities 

Next we compare the flow of various susceptibilities. First 
we analyze the real parts of the static d-wave pairing 
susceptibility Xd-wavc and the AF susceptibility Xs{k,tt). 
Moreover, as a kind of cross-check, we analyze the cou- 
pling to uniform static charge and spin holds for electrons 
at different points on the Fermi surface . 

In the saddle point regime of the hole-doped case the 
overall behavior of these probes together with general fea- 
tures of the RG flow of the couplings indicated the ten- 
dency towards formation of an insulating spin liquid: both 
Xd-wzvc and Xs (i", tt) diverge together in a similar way (see 
Fig. H) and in view of the mutual reinforcement between 
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Fig. 10. Flow of coupling functions (left plots) and susceptibil- 
ities (right plots) in the crossover regime of the electron-doped 
case, t' = — 0.3t and T — 0. The upper plots correspond to 
H = ~0.lt where the FS intersects the US and the AF processes 
processes diverge at non-zero energy scale. The lower plots 
are for n — 0.04i, which is further away from half-filling. The 
dashed lines in the left plots show the flow of Cooper couplings 
at the FS, the solid lines umklapp processes in the BZ diagonal 
Va(4,4,21), Va(3,4,21), and Va(5,4,21). The dashed-dotted 
lines denote Cooper processes between the saddle point regions 
with band energies below the FS (discussed in the text). In the 
right plots the solid line denotes the d-wave pairing susceptibil- 
ity and the dashed-dotted line the AF susceptibility x s (k,tt). 
Note the logarithmic scale axis. The inset in the lower right 
plot shows the Cooper pair scattering V(ki, — ki, k^) close to 
the instability as function of the angle of k3 around the FS. 
The filled squares (open diamonds) correspond to scattering 
between k-space points below (above) the FS at band energy 
±0.8t while the circles connected by straight solid lines show 
the values at for pairs on the FS. ki is at angle w — it below 
(above) or on the FS (see Fig. ^) . 



both channels it is very plausible that the strong coupling 
state will incorporate both types of fluctuations. On the 
other hand the coupling of the saddle point regions to 
uniform charge and spin fields becomes increasingly sup- 
pressed in the RG flow (the charge couplings are shown in 
the left plot in Fig. [ll]) , suggesting the opening of charge 
and spin gaps around the saddle points. In contrast to that 
the charge coupling for quasiparticles in the BZ diagonals 
does not renormalize to zero. This anisotropic flow is con- 
sistent with a truncation of the FS at the saddle points 
while the BZ diagonals remain metallic. The flow of the 
scattering rates described in Sec. ^ corroborates this sce- 
nario. 

In the electron-doped case the situation is rather differ- 
ent. Since d-wave and AF channel do not couple strongly 
the corresponding susceptibilities have very different flows. 
While x s (7r, 7r) already grows considerably at higher scales 
and eventually diverges if the system is not far away from 
half filling, 

Xd-wavc becomes large only at low scales. We 




Fig. 11. Flow of the coupling of quasiparticles at different 
points on the FS to uniform external charge fields, normalized 
to the values with the bare interactions U — 3t. The left plot 
corresponds to the saddle-point regime of the hole-doped case 
fi — —t, T — 0.04£, while the right plot shows data for the 
electron-doped case /i = —O.lt, T — O.Olt. The solid lines are 
for k-points on the FS in the vicinity of the saddle points, the 
dashed lines for k-points on the FS in the BZ diagonal. Af reeze 
is the energy scale down to which the flow of the coupling func- 
tion is taken into account. Below Af reexe the coupling functions 
are kept constant. The insets show the angular variation along 
the FS for different values of A{ reeze . 

therefore expect that the system will have only two dif- 
ferent phases if we allow for a sufficient amount of 3D 
coupling: one AF phase with higher critical scale and tem- 
perature, and, if we increase the band filling such that the 
instability in the latter channel becomes cut off by the FS 
curvature, a rapid drop in the energy scale with a low- 
T Kohn-Luttinger-type d-wave phase. This - compared 
to the saddle point regime - more conventional picture 
is also supported by the behavior of the uniform suscep- 
tibilities and the fc-space dependence of the coupling to 
the corresponding external fields: the coupling to uniform 
charge fields (see right plot in Fig. [ll]) is not strongly 
suppressed, i.e. there are no pronounced indications for 
incompressible or truncated parts of the FS in our weak- 
coupling analysis. The reduction of the charge couplings 
becomes even weaker when the filling is increased. In con- 
trast to the hole-doped case the strongest suppression now 
occurs for wavevectors in the BZ diagonals, in agreement 
with the fact that Umklapp processes flow most strongly 
there (see inset in the right plot of Fig. [ll]). Again this 
type of anisotropy compares well the angular-dependent 
quasiparticle lifetimes obtained with the RG calculations 
described below and FLEX (fluctuation exchange approx- 
imation) 0. 



5 Tentative phase diagram for the 
electron-doped case 

The tentative phase diagram for the hole-doped Hubbard 
model extracted from the one-loop flow to strong coupling 
was discussed in (h]]. Here present a schematic phase di- 
agram for the electron-doped side. 

Generally with this type of calculations it is difficult 
to determine exact phase boundaries, e.g. the parameters 
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where the ground state changes from AF to d-wave super- 
conducting long range order. Nor can we actually prove 
the existence of the phases suggested by one-loop flow. 
Selfenergy and higher order effects will become large once 
the coupling functions exceed the order of the band width 
and from the flow below that scale it is impossible to tell 
whether the different classes of couplings really diverge or 
not. Similarly fluctuations around possible ordered states 
are not taken into account properly. Nevertheless since 
in the electron-doped case we observe two different ten- 
dencies which are weakly coupled and dominate on op- 
posite sides of the considered density range, we can well 
discriminate between distinct regimes with different domi- 
nant fluctuations. If long range order is possible (at T = 
or at finite T if we add some degree of three-dimensional 
coupling), its type will most probably correspond to the 
dominant fluctuations visible in the one-loop RG. Bearing 
in mind these remarks, we present a schematic phase dia- 
gram for the case t' = — 0.3t and initial U = 3t in Fig. [l2| . 
For chemical potential /i < 0.03t (which corresponds to 
(n) w 1.21 per site) we find an AF dominated instability 
with comparably high critical temperature T c (above T c 
we can integrate down to A — with all couplings stay- 
ing smaller than 15£.) . For larger chemical potential the 
instability is d-wave dominated. The criterion we used for 
the distinction was the derivative of Xd-wavc and x s (f, tt) 
with respect to A when the maximal coupling grows larger 
than 12t. We repeat that this distinction does not allow 
any quantitative predictions for the true phase diagram. 
In particular the critical scales and temperatures for the 
d-wave phase strongly depend on the criterion used to sep- 
arate the regimes]]. 

In the electron-doped case considered here there are 
two disparate energy scales. The first one, which can be 
high close to half filling, is related to the AF-(7r, it) pro- 
cesses between the BZ diagonals. The second energy scale 
is given by the critical scale of the d-wave Cooper channel 
and remains much lower close to half-filling. Only if the AF 
energy scale becomes too low at increased doping, when 
the FS curvature prevents a real flow to strong coupling 
of the AF processes, an instability in the d-wave Cooper 
channel can occur. In the present case, the crossover en- 
ergy scale is mainly determined by the particle density 
which regulates the overlap of the FS with the US, i.e. the 
low energy phase space available for elastic scattering with 
momentum transfer {%, w). A qualitatively similar picture 
was found by Manske et al. [19] who solved the general- 
ized Eliashberg equations with a spin-fluctuation induced 
pairing interaction obtained within a FLEX scheme. 

As mentioned above, on the electron-doped side the 
flow to strong coupling is qualitatively similar for all t' < 
0, only the critical scales and relative widths of the regimes 
change. Thus for band fillings larger than one particle per 



1 If we chose to compare the susceptibilities at higher values 
of the couplings, i.e. 20t instead of 12t, the d-wave phase would 
extend further to smaller electron-dopings (to approx. /i = 
O.Olt). Very close to the instability however, when the coupling 
function becomes huge, our one-loop RG scheme cannot be 
justified. 
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Fig. 12. Dependence of the flow to strong coupling on /j, and 
T in the electron-doped regime for t' — — 0.3t and U = 3t. 
The thick line denotes the critical temperature above which 
we can integrate down to A = with all couplings remain- 
ing smaller than 15t. In the d-wave dominated regime Xd-wave 
grows stronger, while in the AF regime Xs( 7r , i") dominates (the 
precise criteria are given in the text) . For the displayed values 
of n the electron-doping varies between 0.2 (fi — — O.lt) and 
0.23 (At = 0.08*). 



site our observations for t' = — 0.3t are in full qualitative 
agreement with the results of Zanchi and Schulz Q for 
t' = and Halboth and Metzner |||| for smaller absolute 
values of the next-nearest neighbor hopping t' . In particu- 
lar in the analysis of Ref. || the d-wave channel becomes 
dominant roughly when the FS loses its intersection with 
the US upon increasing the electron density. 



6 Scattering rates 
6.1 Hole-doped case 

Here we describe our results for the quasiparticle scatter- 
ing rates in the hole-doped case. There we are mainly in- 
terested in electron densities corresponding to the saddle 
point regime |l0| briefly described above, because there 
the flow to strong coupling becomes qualitatively differ- 
ent from the d-wave superconducting instability at lower 
densities. Again, since we cannot prove the FS truncation 
suggested by the flow in the saddle point regime, our main 
goal is to analyze whether the behavior of the scattering 
rates is consistent with this scenario. 

Typical results for the scattering rates in the saddle 
point regime (/i = — Lit) above the critical temp era ture 



obtained with the method described in Subsec. 2.2 are 
shown in Fig. |l3|. The temperature was chosen such that 
the couplings Va do not become too large when we in- 
tegrate the flow to A = 0, for the lowest T shown the 
maximum coupling reaches w 11.5t. This is already larger 
than the bandwidth, but this high value is only reached at 
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low scales and most couplings remain smaller. Therefore 
the results should be qualitatively correct. An anisotropy 
in the scattering rate is clearly observable but not too pro- 
nounced, the maximum ratio between Im S at the saddle 
points and in the diagonals at 7r/4 is rs 2. The tempera- 
ture dependence of Im S{kp,ijj — 0) is shown in Fig. [l3] 
for fi = — l.lf which is above the saddle point regime and 
close to the van Hove doping at /i = — 1.2t. All curves 
for the different positions on the Fermi surface show an 
almost linear increase with T. Both, the angular variation 
and the temperature dependence, are qualitatively con- 
sistent with the ARPES results by Valla et al. fljQ. In 
comparison with the model assumptions in fl3|| we note 
that our calculations at temperatures above the flow to 
strong coupling yield a smaller anisotropy in the scatter- 
ing rates and no T 2 behavior for the quasiparticles in the 
BZ diagonal. Our data resemble more FLEX results ob- 
tained by Kontani et al. Jl7| Altmann et al. fl^| . The 
first group was able to describe resistivity and Hall co- 
efficients in qualitative agreement with the experiments. 
For high temperatures the qualitative similarity between 
FLEX calculations and our RG analysis is not acciden- 
tal because for T > T c when the flow does not diverge, 
the main scattering occurs in the (n, tt) spin fluctuation 
channel which the typical FLEX schemes focus on. 

From a theoretical point of view the linear- T depen- 
dence is not unexpected because already the bare two-loop 
diagram with unrenormalized couplings yields this behav- 
ior if the band filling is sufficiently close to the van Hove 
filling Hence it should be considered as an effect 

of the van Hove singularities and does not immediately 
imply a breakdown of the quasiparticle concept . 

For lower temperatures the coupling function flows to 
strong coupling at a non-zero critical scale A c . Therefore 
we cannot integrate the flow down to zero scale and do 
not obtain a good approximation for Im S with the above 
method. At A c the selfenergy diverges together with the 
couplings and we expect that the quasiparticles will be at 
least partially destroyed. Although we cannot really ac- 
cess this strong coupling region, we can gain some insight 
on how this quasiparticle destruction takes place by con- 
sidering the change in the two-loop selfenergy through the 
flow of the coupling function. The results for \x — —t and 
T = 0.04i are shown in Fig. 14. If we freeze the flow of the 



four point vertices Va already at high scales, we basically 
obtain the bare two-loop selfenergy when we integrate the 
RG equation for the selfenergy down to A = 0. There, 
the anisotropy between saddle point regions and BZ di- 
agonals is small. If we now subsequently include the flow 
of Va by reducing ylf reeze , the scattering rate for quasi- 
particles around the saddle points grows strongly. For the 
quasiparticles in the BZ diagonals the scattering rate is 
not much affected by the flow of the couplings. Thus the 
flow of the scattering rates for these values of the band 
filling is consistent with a breaking up of the FS into two 



2 Valla et al. |15| find a leveling off of the scattering rate very 
close to the saddle points at lowest T in the normal state, which 
is not present in our results. If this is an intrinsic effect of the 
Cu02-planes, it may be beyond our weak coupling calculation. 
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Fig. 13. Left: Temperature dependence of imaginary part of 
the selfenergy at four different positions on the FS for the hole- 
doped case with \i — —l.lt and t! = — 0.3i. The dashed lines 
belong to points in the BZ diagonal, the solid lines to points 
closer to the saddle points. Right: Imaginary part of the selfen- 
ergy versus angle around the FS in the saddle point regime at 
jj, = —l.lt for different temperatures above the instability. The 
lower line (circles) corresponds to T — 0.154, the middle line 
(diamonds) to T = 0.184 and the upper line is at T = 0.214. 
For T = 0.154 the largest coupling at A — is « 11.5i, for 
T — 0.214 they reach » 8.94. The saddle points are at angles 
and 7r/2. 
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Fig. 14. Left: Change of the scattering rate through the flow of 
the interactions Va for different positions on the Fermi surface 
in the hole-doped case, calculated on a 48 x 3 system. /lf rcozo 
is the scale below which the interactions are kept constant. 
Right: Angular variation of Im £(\tF,iJ = 0) with Icf varying 
on the FS from one to the next saddle point. The different lines 
shows E(kp,u) = 0) with the couplings stopped at V/i jmax = 64 
(bottom line), 94, 12t and 154 (top line). The inset shows the 
12 points on a quarter of the FS. For both plots fi = —l.lt and 
T — 0.044 corresponding to the saddle point regime. 



distinct regions, as suggested in |l0| based on an analysis 
of the flow of the coupling function and susceptibilities: 
around the saddle points, the quasiparticles are subject 
to divergent scattering processes, there we also found the 
strong suppression of the charge compressibility. In the 
BZ diagonals which appeared to remain compressible, the 
scattering rate stays in the weak coupling range, consis- 
tent with the FS remaining untruncated in these regions. 

When the band filling is increased further towards half- 
filling the energy and temperature scale for the flow to 
strong coupling rises and the anisotropy of the scattering 
rates at temperatures above this flow to strong coupling 
becomes weaker. The reason is that the scattering pro- 
cesses between the saddle point regions lose their dom- 
inant role, and the strongest scattering processes occur 
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between FS parts connected by (tt, tt), now closer to the 
BZ diagonals. 



6.2 Electron-doped case 

In Fig. [l5| we show the results for the electron-doped 
case obtained at higher temperatures where the flow does 
not exceed the order of the bandwidth. There are strik- 
ing differences compared to the results obtained in the 
hole-doped saddle point regime. First the overall magni- 
tude of Im E{\s.p,u) = 0) at a given temperature is much 
smaller, almost by an order of magnitude. We attribute 
this to the reduced density of states around the FS, the 
van Hove singularities are far away from the FS. Second 
the temperature dependence of the scattering rates is far 
from being linear and appears to be roughly consistent 
with a T 2 -behavior at low temperatures above the flow 
to strong coupling. The data for k-points in the BZ di- 
agonal do not extrapolate to for T — > due to the in- 
crease of the coupling function at lower T. Further, the 
anisotropy of the scattering rate is different and much 
weaker than in the hole-doped case. The FS parts which 
feel the strongest scattering at lower T are now in the BZ 
diagonal (see right plot in Fig. |l5| for angular dependence). 
This is consistent with the flow of the interactions which 
grow fastest in these regions. Our results also agree quali- 
tatively with FLEX calculations p7|,^9[, which show that 
the in the electron-doped case the "hot spots" with the 
shortest quasiparticle-lifetime move towards the BZ diag- 
onal. Nonetheless due to the weakness of the anisotropy 
the data do not suggest a truncation of the FS like in the 
hole-doped case. This can be also seen at lower tempera- 
tures (see Fig. 16]), where flow of the scattering rates re- 
mains comparably weak when the coupling function flows 
to strong coupling. 



7 Discussion 

We have presented iV-patch RG results for the two-dimen- 
sional t-t' Hubbard model with particle density n > 1 
per site and t' = — 0.3t. We have shown that the RG 
flow to strong coupling is qualitatively different from the 
n < 1 case because of the different location of the FS in 
both cases. The saddle-point regime jl(| which we asso- 
ciated with the pseudogap state in the (hole-)underdoped 
cuprates is absent in our RG study of the electron-doped 
side. Instead the picture closely resembles the observations 
made by Zanchi and Schulz || and Halboth and Metzner 
|^,D for zero or smaller absolute value of the next-nearest 
neighbor hopping t' . Upon increasing the electron den- 
sity we find a rapid decrease of the critical scale for the 
instability where the flow of the otherwise strongly grow- 
ing AF susceptibility is cut off. This leaves a small win- 
dow for a Kohn-Luttinger-type d x 2_ y 2-w&ve instability at 

3 In an isotropic two-dimensional Fermi liquid, we expect 

im E(u = o,k F ) ~ (r/r F ) 2 iog(r F /r) §§,§if. Since our cal - 

culations are rather qualitative we do not attempt a detailed 
comparison of the data with this T-dependence. 
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Fig. 15. Left: Temperature dependence of imaginary part of 
the selfenergy at four different positions on the FS for the 
electron-doped case. The dashed lines belong to points in the 
BZ diagonal, the solid lines to points closer to the saddle points. 
Right: Angular variation of Im S(k-,uj = 0) with k F varying 
on the FS from one to the next saddle point. The different lines 
show E(k.-,u> = 0) at different temperatures from T — O.ldt 
(bottom line) to T = 0.28t (top line). The results in both plots 
were obtained using a 32 x 3 system with t' — — 0.3t and fi = 0. 
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Fig. 16. Left: Change of the scattering rate through the flow of 
the interactions Va for different positions on the Fermi surface 
in the electron-doped case. Af ree ze is the scale below which the 
interactions are kept constant. Right: Angular variations of Im 
S(kp, u) — 0) with k F varying on the FS from one to the next 
saddle point. The different lines shows E(kw,u = 0) with the 
couplings stopped at F/i.max = 6t (bottom line), 9t, 12t and 
15t (top line). The inset shows the 12 points on a quarter of 
the FS. For both plots fi = ~0.lt and T = 0.04i. 



lower scales and temperatures. The tendency towards in- 
compressibility and truncation of parts of the FS is weaker 
than in the hole-doped case and strongest in the BZ diag- 
onal. At temperatures above the strongly coupled phase 
we find slightly anisotropic scattering rates with strongest 
scattering of quasiparticles, i.e. "hot spots" in the BZ diag- 
onal. The temperature dependence of the scattering rates 
is roughly quadratic. 

Indeed the experimental phase diagram of the electron- 
doped cuprates shows a relatively wide AF phase with 
a boundary to superconducting phase with lower critical 
temperature. There is increasing experimental evidence 
^4|[2^,^6|,^7| that the pairing symmetry of this phase is 
d x 2_ y 2. Moreover the in-plane resistivity in the normal 



29 



state shows a quadratic temperature dependence 

Note that in our RG study the difference between the 
energy scales characteristic to these phases is considerably 
larger than the experimental values. However one should 
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keep in mind that our one-loop analysis can only provide a 
qualitative picture for the understanding of the real mate- 
rials. In particular we cannot make any predictions about 
the stability of the phases suggested by the RG flow and 
the precise location of the transition from AF to d-wave 
regime. In fact the Fermi surface in Ndi.ssCeo.isCuO^ 
with a superconducting T c of 2AK seen by ARPES Q 
still intersects the Umklapp surface. Thus it might be that 
our criterion for the transition from AF to d-wave regime 
overestimates the width of the AF regime. 

The saddle point regime in the hole-doped case with 
the diverging umklapp processes was interpreted as a pre- 
cursor of the Mott state where the FS becomes truncated 
at the saddle points . The anisotropic scattering rates 
described above, which grow strongest at the saddle points 
as the flow goes to strong coupling, are consistent with 
this scenario. Due to the location of the FS the diverg- 
ing umklapp processes also generate strong d-wave pairing 
correlations. In this situation the d- wave-paired conden- 
sate which may form at sufficiently low temperature will 
still be subject to the umklapp scattering at the saddle 
points where the FS becomes truncated. This may result 
in a reduced superfluid weight arising from the open parts 
of the FS in the BZ diagonals only. 

In the electron-doped case the umklapp and other AF 
processes predominantly occur between the BZ diagonals 
and do not strongly drive the <i-wave correlations, they 
can only lead to a rapid growth of AF correlations. Thus, 
instead of being strongly coupled and mutually reinforc- 
ing each other in the flow, AF and d-wave channel have 
disparate characteristic energy scales and influence each 
other only weakly. As soon as the umklapp processes in 
the BZ diagonals flow to strong coupling, the instability 
is entirely dominated by AF correlations because the en- 
ergy scale of the d-wetve instability is much lower. This 
implies that coming from the electron-overdoped side the 
precursors of the Mott state only occur in the AF phase 
and should have less influence on the d-wave state. In par- 
ticular we speculate that in the superconducting state all 
electrons contribute to the superfluid weight in contrast 
with the hole-(under)doped case. 



It is a pleasure to thank T.M. Rice, M. Salmhofer, U. 
Ledermann, K. Lehur, M. Zhitomirsky, R. Hlubina and D. 
Vollhardt for stimulating discussions. The Swiss National 
Science Foundation is acknowledged for financial support. 
The numerical calculations were performed on the Asgard 
Beowulf cluster of ETH Zurich. 



10. 

11. 
12. 
13. 
14. 
15. 

16. 



17. 



18. 



20. 
21. 

22. 
23. 

24. 

25. 

26. 
27. 



28. 



29. 



P. Lederer, G. Montambaux and D. Poilblanc, J.Phys. 48, 
1613 (1987). 

D. Zanchi and H.J. Schulz, Europhys.Lett. 44, 235 (1997); 
D. Zanchi and H.J. Schulz, Phys. Rev. B 61 13609 (2000). 
N. Furukawa, T.M. Rice, and M. Salmhofer, 
Phys.Rev.Lett. 81, 3195 (1998). 

C.J. Halboth and W. Metzner, Phys.Rev.B 61, 7364 
(2000). 

C.J. Halboth and W. Metzner, Phys. Rev. Lett. 85, 5162 
(2000). 

C. Honerkamp, M. Salmhofer, N. Furukawa, and T.M. Rice, 
Phys. Rev. B 63 35109 (2001). 

Shan- Wen Tsai and J. B. Marston , |cond-mat/001030C . 
M.P.A. Fisher, :ond- mat/9806164 and references therein. 
L.B. Ioffe, A.J. Millis, Phys. Rev. B 58, 11631 (1998). 
R. Hlubina, Phys. Rev. B58, 8240 (1998). 
T.Valla, A.V.Federov, P.D.Johnson, Q.Li, G.D.Gu and 
N.Koshizuka, Phys. Rev. Lett, 85, 828 (2000). 
M. Salmhofer and C. Honerkamp, to be published in Prog. 
Theo. Physics 105 (2001); for a related formalism, see also 
M. Salmhofer, Renormalization: An Introduction, Texts 
and Monographs in Physics, Springer (1999). 
H.Kontani, K.Kanki, and K Ueda, Phys.Rev.B 59, 14723 
(1999). See also H. Kontani, |cond-mat/0011327 . 

D. Manske, I. Eremin, and K.H. Bennemann, Phys Rev. 
62, 13922 (2000). 



19. J. Altmann , W. Brenig and A. P. Kampf, cond- 



mat/9707267 (unpublished). 



PA. Lee and N. Read, Phys. Rev. Lett. 58, 2691 (1987). 
G.Kastrinakis, Physica C 318, 497 (2000); Physica C 340, 
119 (2000) . 

A.V. Chaplik, Sov. Phys. JETP 33, 997 (1971). 

C. Hodges, H. Smith, and J.W. Wilkins, Phys.Rev. B 4, 

302 (1971). 

C.C. Tsuei, JR. Kirtley, Phys. Rev. Lett. 85, 182 (2000); 
Rev. Mod. Phys. 72, 969 (2000). 

R. Prozorov, R. W. Giannetta, P. Fournier, and R. L. 

Greene, Phys. Rev. Lett. 85, 3700 (2000). 

A. Mourachkine, Europhys. Lett. 50, 663 (2000). 

N.P. Armitage, D.H. Lu, D.L. Feng, C.Kim, ADamascelli, 

K.M. Shen, F. Ronning, Y. Onose, Y. Taguchi, Y. Tokura, 

and Z.-X. Shen, Phys.Rev.Lett. 86, 1126 (2001). 

C.C. Tsuei, A. Gupta, and G. Koren, Physica C 161, 415 

(1989). 

J.L. Peng,E. Maiser, T. Venkatesan, R.L. Greene, and G. 
Czjzek, Phys.Rev. B 55, R6145 (1997). 



References 

1. P.W. Anderson, The Theory of Superconductivity in the 
High-T c Cuprates, Princeton University Press (1997). 

2. For a review see: D.J. Scalapino, J. Low Temp. Phys. 117, 
179 (1999). 

3. I. Dzyaloshinskii, Sov.Phys.JETP 66, 848 (1987); I. 
Dzyaloshinskii, J.Phys.I (France) 6, 119 (1996). 

4. H.J. Schulz, Europhys.Lett. 4, 609 (1987). 



